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THE  IMPLICIT  COMPLEMENTARITY  PROBLEM:  PART  II 
Jong-Shl  Pang 

ABSTRACT.  In  Part  I  of  this  study,  we  have  defined  the  implicit  - p1 - 

tarlty  problem  and  investigated  its  existence  and  uniqueness  of  solution. 

In  the  present  paper,  we  establish  a  convergence  theory  fer  a  certain 
iterative  algorithm  to  solve  the  implicit  complementarity  problem.  He  also 
demonstrate  how  the  algorithm  includes  as  special  cases  many  existing 
iterative  methods  for  solving  a  linear  complementarity  problem. 


Key  Words.  Implicit  complementarity 


Iterative  method 


Successive  overrelaxation 
Convergence  theory 


1*  Introduction.  Given  an  n  by  n  matrix  A,  n-vector  b  and  a  mapping  m  from 
Ra  Into  Itself,  the  Implicit  complementarity  problem  (ICP) ,  denoted  by  the 
triple  (A,b,m)  la  to  find  a  vector  x  In  Rn  satisfying 

Ax  +  b  £  0,  x  >  m(x)  and  (Ax  +  b)T(x  -  m(x))  -  0. 

In  a  previous  paper  [11],  we  have  shown  how  various  complementarity  problems 
can  be  cast  as  an  ICR  and  studied  the  existence  and  uniqueness  of  solution 
to  the  ICP.  The  fundamental  tool  employed  In  our  study  is  a  certain 
implicitly  defined  mapping  F.  With  this  mapping  F,  we  generate  a  sequence 
of  vectors  {u  }  iteratively  by 

(1)  -  Ftu1*)  k  >  0 

where  u°  la  a  given  initial  vector. 

In  this  paper,  we  study  the  convergence  of  the  iterative  scheme  (1) 
to  a  solution  of  the  ICP.  Our  purpose  is  twofold.  First,  we  show  that 
many  iterative  methods  for  solving  the  linear  complementarity  problem  can  be 
unified  and  extended  under  this  general  scheme  (1).  Second  we  establish 
a  theory  for  the  convergence  of  the  scheme. 

Over  the  past  several  years,  there  has  been  an  increasing  amount  of 
studies  on  Iterative  methods  for  solving  the  linear  complementarity  problem. 
Several  recent  references  are  [1,  2,  3,  4,  5,  7].  In  all  except  two  [1,  2) 
of  these  previous  works,  the  convergence  proofs  of  the  methods  rely  heavily 
on  the  synmetry  of  matrix  involved.  Part  of  the  contribution  of  the 
present  research  is  that  we  have  provided  a  general  framework  for  the  study 
of  these  iterative  methods  for  the  linear  complementarity  problem.  More 
Importantly,  our  method  of  convergence  proof  is  based  on  a  theory  of 
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contraction  mappings  and  does  not  rely  on  matrix  synmetry.  (This  same 
approach  vas  used  by  Aganaglc  and  Ahn  In  their  studies.) 

The  rest  o£  the  paper  Is  organised  in  two  sections.  In  the  next 
section,  after  briefly  reviewing  some  background  materials  (Including  the 
definition  of  the  mapping  F) ,  we  show  how  various  iterative  methods  for  the 
linear  complementarity  problem  can  be  obtained  as  special  cases  of  the  general 
scheme  (1).  In  the  last  section,  we  develop  a  convergence  theory  for  the 
iterative  scheme  (1)  and  discuss  some  of  its  specializations. 

2.  Iterative  Methods  for  the  LCP.  Given  an  n  by  n  matrix  A  and  n-vector 
b,  the  linear  complementarity  problem  (LCP),  denoted  by  the  pair  (A,b)  is 
a  special  case  of  the  ICP(A,b,m)  where  m  is  the  zero  mapping. 

The  mapping  F  used  to  define  the  iterative  scheme  (1)  is  constructed 
indie  following  way.  See  [11].  Given  a  P-splittlng  (B,  C)  of  the  matrix 
A,  l.e.  A  ■  B  -  C  where  B  Is  a  P -matrix  (l.e.  has  positive  principal  minors), 
for  each  vector  u  €  R°,  F(u)  is  the  unique  solution  to  the  LCP 

Bx  +  (b  -  Cu)  >0,  x  >  m(u),  and  (Bx  +  b  -  Cu)T(x  -  m(u))  ■  0. 

(This  latter  problem  Is  not  an  LCP  In  the  ordinary  sense.  However,  by  the 
obvious  translation  of  variables  y  •  x  -  m(u),  it  becomes  the  LCP(B,  b  -  Cu 
+  Bm(u))).  It  is  easy  to  see  that  a  vector  u*  Is  a  fixed  point  of  F  If  and 
only  if  It  Is  a  solution  of  the  IC?(A,b,m).  For  more  discussion  on  this 
mapping  F  and  its  role  in  the  study  of  the  ICP,  see  [11].  (The  P-splittlng 
used  to  define  the  mapping  F  should  not  be  confused  with  the  P-regular 
splitting  of  a  matrix  used  frequently  In  the  numerical  analysis  literature 
(see  Ortega  [8]  e.g.).  For  more  discussion  on  various  matrix  splittings 
smployed  in  the  study  of  the  ICP,  see  [11]). 


Many  iterative  methods  for  solving  the  LCP  are  of  Che  successive 
overrelaxation  (SOR)  type.  These  methods  are  based  on  their  counterparts 
for  solving  systems  of  linear  equations.  The  following  algorithm,  due  to 
Cryer  [5]  is  the  modified  point  SOR  method  for  solving  the  LCP  (M,  q) 
where  the  matrix  M  is  n  by  n. 

Algorithm  I.  (Cryer)  Let  z®  be  an  arbitrary  nonnegative  vector  and  id*  a 
scalar  in  the  interval  (0,  2).  Generate  the  sequence  fzK}  as  follows. 

For  i  -  l,...,n,  let 

M-r + j  >  i 

and 

k+1  f  k  .  .  . «k+l  k.  t 

zt  -  max[0,  zi  +  m*< z±  -  zp}. 

In  f4j,  Cottle,  Golub  and  Sacher  extended  the  above  point  method  to  a 
block  iterative  scheme.  More  precisely,  let  the  matrix  M  be  partitioned 
into  submatrices  (i,j  ■  l,...,n)  where  is  of  order  n^  by  n^.  The 
vectors  z  and  q  are  partitoned  accordingly.  The  algorithm  below  is  the 
modified  block  SOR  method  for  solving  the  LCP  (M,  q) .  The  point  version  of 
the  algorithm  corresponds  to  the  case  where  each  block  size  n^  is  equal  to 
one. 

Algorithm  II.  (Cottle,  Golub  and  Sacher)  Let  z°  -  (z®,...,z®)  be  an 

arbitrary  nonnegative  vector  and  uj*  a  scalar  in  the  interval  (0,  2) . 

k  k  u 

Generate  the  sequence  £z  “  (z£,  ...,zp)  as  follows.  For  1  -  l,...,n, 

let  £^+1  solve  the  LCP  (M. . ,  z  +  £  M  z^1  +  £  M.  .zS  and  let 

J  <  i  J  J  j  >  i  3  3 
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s 


fe+1 

1 


k  ,  lc+1 .  «k+l 
*i  +  u,l  <*i 


where 


<»k+l  -  max{(i)  :  u>  <  va*  ,  z£  +  u»(*^+l  -  z£)  >  0). 

It  was  pointed  out  in  [4]  that  if  w*  €(0,  1],  then  ui^1  -  ta*  for 
each  1  and  each  kj  on  the  other  hand,  if  to*  €(1,  2),  then  wk+1  >  l. 

Moreover,  if  a*  €(1,  2),  then  Algorithm  I  and  the  point  version  of 
Algorithm  II  are  equivalent  in  the  sense  that  the  same  sequence  of 
iterates  (z  ]  is  generated. 

Motivated  by  several  earlier  works,  including  that  of  Cryer,  Mangas arisn 
[7]  proposes  the  following  fairly  general  iterative  algorithm  for  solving 
the  LCP  (M,  q) . 


Algorithm  III.  (Mangas arian)  Let  z*  be  an  arbitrary  nonnegative  vector, 
\€(0,  1)  and  to*  >  0.  Generate  the  sequence  [zkJ  as  follows.  Let 

zk+1  “  A  max  Co,  zk  -  ui*E(q  +  Mzk  +  KCz1'*1  -  a1*))}  +  (1  -  \)zk 

where  E  is  a  poaltive  definite  diagonal  matrix. 

In  its  original  formulation,  the  matrices  E  and  K  are  allowed  to 
v«y  Itom  one  Iteration  to  the  next.  For  our  purpose  here,  they  are 
chosen  to  be  fixed  throughout  the  algorithm.  In  the  reference. 

Mangas arian  commented  that  in  order  for  the  algorithm  to  be  pratlcal, 
the  matrix  K  must  be  either  strictly  upper  or  lower  triangular.  As  we 
shall  show  later,  this  restriction  on  the  choice  of  K  can  be  relaxed 
considerably. 
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Recognislng  that  the  LCP  can  be  formulated  as  e  certain  fixed  point 
problem,  Aganaglc  [1]  proposed  the  following  iterative  algorltim  for 
aolving  the  LCP  (M,  q) .  Although  for  the  convergence  of  this  algorithm, 
the  initial  vector  c  need  not  be  chosen  nonnegative j  in  order  to  compare 
with  the  above  three  algorithms,  we  restrict  Aganaglc' s  algorithm  to 
start  with  a  nonnegative  vector. 

Algorithm  IV.  (Aganaglc)  Let  s®  be  an  arbitrary  nonnegative  vector, 

\€[0,  1]  and  w*  >  0.  Generate  the  sequence  fa*}  as  follows.  Let 

a**1  -  xk  -  X  min{u»*(MSk  +  q),  sk]. 

It  is  easy  to  see  that  both  Algorithms  I  and  IV  are  special  cases  of 
III.  In  fact,  with  X  ■  1,  E  ■  D  and  K  *  L  where  D  and  L  are  respectively 
the  diagonal  end  strictly  lower  triangular  parts  of  the  matrix  M, 

Algorithm  III  reduces  to  I;  whereas  with  K  ■  0  and  E  •  I,  Algorithm  III 
becomes  IV.  (Rigorously  speaking,  the  diagonal  entries  of  the  matrix  M 
need  to  be  positive  in  order  for  Algorithm  III  to  include  I  as  a  special 
case.  In  Cryer's  original  paper  [5],  this  requirement  is  always  met 
because  M  is  assumed  to  be  symmetric  positive  definite.)  In  order  to  see 
how  Algorltham  II  and  III  are  related,  we  establish  the  next  two  prop¬ 
ositions  which  describe  an  equivalent  way  for  generating  the  sequence  of 
iterates . 

Proposition  2.1.  Suppose  that  each  diagonal  block  Mlt  of  the  matrix  M  is 
a  P-satrlx.  Let  D,  L  and  U  be  respectively  die  block  diagonal,  the  block 
strictly  lower  and  strictly  upper  triangular  parts  of  M  respectively.  Then 
provided  that  either  each  n^  is  equal  to  1,  or  the  parameter  os*  is  in 
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(0,  1],  the  sequence  {*  ' j  generated  by  Algorithm  II  can  be  obtained  In  the 

k  M 

following  way.  Given  a  nonnegatlve,  let  z  be  the  unique  solution  to  the 
LCP 

(2)  u  -  <tf*q  +  [(0*0  -  (1  -  (tfAODjz*6  +  (D  +  w*L)z  >  0 
v  -  -  max(Ot  1  -  <o*)z'C  +  z  >  0,  u*v  •  0. 


Proof.  We  shall  prove  the  theorem  only  for  the  case  where  each  n^  is 
equal  to  1.  The  other  case  can  be  proved  in  a  similar  way.  First  of  all, 
observe  that  if  each  n^  is  equal  to  1,  then  for  each  1,  z£+*  »  maxfo,  z£ 

+  ob*(z£+1  -  z£)}.  Moreover, 


*t  +  )  *  max  {(1  -  <»*)*£, 

(1  -  <o*)zjj  -  «o*(qi  +  jEi  m^zj**1  + 
It  therefore  follows  that 


k+1 


fmax(0,  (1  -  (0*)zj  )  , 


_k+l 


(1  -  (0*)z  -  co*(q  +  E  m  z  +  E  m..z  ) 

1  1  j  <  i  lJ  J  j  >  i  1 J  J  11 

or  equivalently, 

u.  *  <o*q.  +  [u*  2  -  (1  -  (0*)m  . z*]  +  (m  .z^1  +  otf*  E  m. .a**1)  >  0 

1  1  j  >  i  U  J  ii  i  11  1  j  <  i  J  • 


k  Uff 

v^  »  -  max[0,  (1  -  (0*)z^]  +  z^  >0  and  ui  *  vi  " 

k  k  k 

Since  z  is  nonnegative,  aax[0,  (1  -  to*)z  ]  •  [max(0,  1  -  <o*)]z  .  Hence  the 

vector  z  is  a  solution  to  the  LCP (2).  The  reason  that  it  is  the  unique 

solution  is  because  the  aatrlx  D  +  oj*L  is  a  P -matrix  by  the  assumptions  on 


D  and  L 


Q.B.D 
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Ui.1 

Propoiition  2.2.  The  vector  z  generated  by  Algorithm  III  Is  e  solution  to 
the  LCP 

(3)  u  “  Xui*q  +  (XtU*(M  -  K)  -  D)*k  +  (D  +  XuJ*K)z  >  0 
v  -  -  (1  -  X)zk  +  z  >  0  ,  uTv  «  0, 
where  D  “  E~^. 

Proof.  By  definition, 

zk+1  -  max  [(1  -  X)zk,  zk  -  X<u*E(q  +  (M  -  K)zk  +  Kzk+1) }. 

fcfi 

Thus,  z  is  a  solution  to  the  problem 

u  -  Xuj*Eq  +  (Xx»*E(M  -  K)  -  I)zk  +  (I  +  Xcu*EK)z  >  0 
v  -  -  (1  -  X)zk  +  z  >  0  ,  uTv  “  0. 

Since  E  is  positive  diagonal,  this  latter  problem  is  obviously  equivalent 

to  <3>-  Q.E.D. 

It  is  important  to  point  out  that  Proposition  2.2  holds  regardless 
of  what  the  matrix  K  is.  Consequently,  a  sufficient  condition  for  the 
sequence  (s  ]  to  be  well-defined  in  Algorithm  III  is  that  (D  +  Xto*K)  is 
a  P-matrix.  This  condition  is  certainly  satisfied  if  K  is  strictly 
triangular  and  if  D  is  a  positive  definite  diagonal  matrix. . 

From  the  two  Propositions  2.1  and  2.2,  it  should  be  obvious  that  the 
point  version  of  Algorithm  II  is  a  special  case  of  Algorithm  III.  In 
fact,  if  the  parameter  os*  (in  Algorithm  II)  is  in  (0,  1],  then  by  setting 
in  Algorithm  III,  X  ■  w*  (of  II),  V*  (of  III)  -  1,  and  K  -  L,  we  obtain 
Algorithm  II.  On  the  other  hand,  if  v*  (of  II)  exceeds  one,  then  by 
setting  la  Algorithm  III,  X  ■  1  and  K  ■  L,  we  obtain  Algorithm  II  as  wall. 


•8» 


For  given  scalars  and  kg  with  ^  >  0  and  <  1,  the  ICP  (k^M,  kgq,  a) 
where  m(z)  *  k^x  for  all  x,  is  obviously  equivalent  to  the  LCP  (M,  q). 
Consider  the  splitting  (B,  C)  of  the  matrix  k^M  : 

(4)  B  *  D  +  k^K  and  C  *  D  +  k^[E  •  M] 

where  O  and  K  are  arbitrary  matrices  such  that  B  is  a  P-matrix.  It  is  easy 
to  see  that  by  appropriately  choosing  the  scalars  k^  and  kg  and  the  matrices 
D  and  K,  the  general  iterative  algorithm  defined  by  (1)  for  the  ICP  reduces 
to  Mangasarisn's  algorithm  III  (which  includes  Algorithms  I  and  IV  and  the 
point  version  of  Algorithm  II)  as  well  as  to  the  general  block  version 
(with  parameter  <u*  not  exceeding  one)  of  Algorithm  II.  It  is  not  clear  to 
the  author  at  this  stage  whether  Algorithm  II  with  oj*  larger  than  one  can 
also  be  obtained  as  a  special  case  of  the  iterative  scheme  (1). 

3.  A  Convergence  Theory.  Given  an  n  by  n  real  matrix  B,  we  define  its 
comparison  matrix  B  by 

®ii  *  lBul  and  Btj  -  -|BtjI  for  i  ^  j  , 

(see  Varga  [13]).  The  matrix  B  is  said  to  be  an  H-matrix  if  its  comparison 
matrix  B  is  a  P-matrix  (see  Ostrowskl  [9]).  We  call  the  splitting  (B,  C) 
of  the  matrix  A  an  H-splitting  if  B  is  an  H-matrix  with  positive  diagonals. 
Since  any  H-matrix  with  positive  diagonal  entries  is  a  P-matrix  (see  Pang 
[10]),  an  H-splitting  is  necessarily  a  P-spllttlng.  Conversely,  a  P- 
spllttlng  (B,  C)  where  B  is  block  triangular  with  each  diagonal  block 
being  an  H-matrlx  with  positive  diagonals,  is  an  H-splittlng.  In 
particular,  the  splitting  (4)  is  an  H-splltting  if  for  instance,  D  is  a 
positive  definite  diagonal  matrix  and  K  is  strictly  triangular. 


In  what  follows,  we  fix  an  arbitrary  H-aplitting  (B,  C)  of  the  matrix 


A  and  study  the  iterative  scheme  defined  by  (1) .  For  any  two  vectors  u 
and  vt  let  e(ut  v)  and  n(u,  v)  be  the  vectors  whose  components  are  the 
differences  jF(u)^  -  F(v)^|  and  jmOO^  -  mCv)^  respectively.  The  result 
below  establishes  a  basic  relationship  between  these  latter  two  vectors. 

It  is  an  improvement  as  well  as  an  extension  of  the  one  derived  by  Ahn  [2] 
for  a  certain  special  case. 

Proposition  3.1.  Let  (B,  C)  be  an  H- splitting  of  the  matrix  A.  Then 
(5)  e(u,  v)  <  B  1  max[5n(u,  v)  ,  |c|  |u  -  v|], 

where  5  is  the  diagonal  matrix  whose  diagonal  entries  are  those  of  B,  |c| 
and  |u  -  v|  are  respectively,  the  matrix  and  vector  whose  components  are 
the  absolute  values  of  those  of  C  and  u  -  v. 


Proof.  Consider  an  arbitrary  index  i  and  suppose  that 
jF(u)1  -  F(v)iJ  -  F(u)t  -  F(v)i. 

If  F(u)t  ■  m(u)^,  then 


I F<u)i  -  F(v)i|  <  m(u)i  -  m(v)  <  (  m(u)  -  m(v>  |  +  S  (fi  |  |F(u).  -  F(v)  l/B 

j?*i  J  3  J 

On  the  other  hand,  if  fBF(u)  +  (b  -  Cu)]^  ■  0,  then  by  letting  to  denote 
the  i-th  row  of  the  matrix  C,  we  have 


|F(u)  -  F(v)  I  <  C-[(b  -  Cu)  +  fi  B  F(u)  ]  +  [(b  -  Cv).  +  £  B.  4F(v)  J}/B,, 

ji*i  13  3  1  iJ  5  ii 

-  [C.(U  -  v)  -  S  B  <F(u)  -  F(v)) ,  J/B 

J*i  ^  11 

<  flc^u  -  v)  |  +  jBtj]  |F(u)j  -  F(v)  jl  }/Btl« 
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Hence,  It  follows  that  In  either  case 

|F<u>±  -  F(v)t|  <  max{Jm(u)i  -  m^oj  ,  Jc^u  -  v)j/Bil} 

l««)j  -  W,|/»u  • 

Consequently,  we  obtain 

B  e(u,  v)  <  max{5n(u,  v)  ,  Jcj  |u  -  vj}  . 

The  desired  Inequality  (5)  now  follows  from  the  fact  that  B  is  nonnegative 
(see  Fiedler  and  Ptak  [6],  e.g.). 

With  the  proposition  above,  we  state  the  following  principal  con¬ 
vergence  result  for  the  iterative  scheme  (1)  to  solve  the  ICP  (A,b,m). 

Theorem  3.2.  Let  (B,  C)  be  an  H-splitting  of  the  matrix  A.  Suppose  that 
there  exists  a  matrix  E  such  that 


(6a) 

n(u,  v)  <  E|u  -  v| 

for  all  u,  v 

_-l 

(6b) 

p (G)  <  1  where  G 

■  B  raax(DE  ,  |c|)  , 

with  p  denoting  the  spectral  radius.  Then  for  any  initial  vector  u°,  the 

sequence  {u  }  generated  by  (1)  converges  to  a  solution  of  the  ICP  (A,b,m). 

Before  proving  the  theorem,  we  point  out  several  remarks.  First, 

the  matrix  E  must  necessarily  be  nonnegative .  Second,  the  initial  vector 
0 

u  is  not  required  to  be  nonnegative.  Third,  the  theorem  may  be  considered 
as  giving  a  set  of  sufficient  conditions  for  the  existence  of  solution  to 
the  ICP. 

Proof  of  Theorem  3.2.  (Similar  to  Ahn  [2])  First  of  all,  observe  that  the 
matrix  G  is  nonnegative.  Conditions  (5)  and  (6a)  imply  that  for  k  £  1, 
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|  fcfl  k|  .  _i  K  K-li 
Ju  -  u  I  <  Gju  -  u  J . 

Since  p(G)  <  1,  it  follows  that 

(7)  lim  \uM  -  u^|  ■  0. 
k  -*  • 

Next,  by  an  inductive  argument,  we  may  deduce 

\uM  -  u°|  <  k  S  1  G^u1  -  u°|  £  (1  -  G)"1^1  -  u°| 
i  -  0 

where  the  last  inequality  follows  from  the  fact  that  the  matrix  G  is  non¬ 
negative  and  p(G)  <1  (see  Ortega  [8,  p.  26]).  Hence  the  sequence  {uk} 
is  bounded  and  thus  has  an  accumulation  point  u*. 

Let  {u^}  be  a  subsequence  converging  to  u*.  Then  (7)  Implies  that 
c  k4+l, 

tu  J  converges  to  u*  as  well.  From  condition  (6a),  it  follows  that  the 
mapping  m  is  continuous.  Hence,  by  passing  the  limit  k^  —  •  in  the 


k  k-1. 


conditions 


ki+1  *  ki+l 


Bu  +  (b  -  Cu  )  >  0  ,  u  >  m(u  l)  ,  (u  -  m(u^))  (Bu^  +  b  -  Cu  )  - 

we  deduce  that  u*  is  a  solution  to  the  ICP  (A,b,m).  Thus  u*  is  a  fixed 
point  of  the  mapping  F.  By  using  Proposition  3.1  again,  we  obtain 

lu**1  -  u*|  <  G|u^  -  u*|  <  ....  <  G^fu*  -  u*|. 

As  p(G)  <  1,  it  follows  that  the  entire  sequence  {u^}  converges  to  u*. 

This  completes  the  proof  of  the  theorem. 

The  argument  used  in  the  above  proof  is  a  contraction-type  reasoning 
to  establish  the  convergence  of  sequences.  Such  argument  is  used  frequently 
in  the  study  of  iterative  methods  for  systems  of  linear  equations  (see 
Ortega  [8]  e.g.). 


Ws 


-12- 


Before  discussing  some  special  cases  of  Theorem  3.2,  we  give  necessary 
and  sufficient  conditions  for  (6a)  and  (6b)  to  hold.  The  next  result  shows 
that  condition  (6a)  is  equivalent  to  the  fact  that  the  mapping  m  is 
Lipschltz  continuous,  i.e.  there  exists  a  scalar  t*  such  that 

(8)  ||m(u)  -  m(v)j[  <  d||u  -  v|f  for  all  u,  v 

where  the  double  vertical  lines  denote  a  certain  norm  of  Ra. 


Proposition  3.3.  Condition  (6a)  is  satisfied  for  some  matrix  E  if  and  only 
if  the  mapping  m  is  Lipschltz  continuous. 


Proof.  Suppose  that  m  is  Lipschltz  continuous.  With  no  loss  of  generality, 

we  may  assume  that  the  norm  in  (8)  is  the  “-norm  ;  llxl)  »  max  lx,,  I  . 

1  S  i  £  n  * 

Then  for  any  u  and  v,  we  have 

i  JT*  Jm(u) i  '  m(v)J  2  |u  -  v  | 

1  S  1  S  n  x  1  j  «  1  1  1 


so  that 


|m(u)  -  m(v)J  <  p  E |u  -  vj 


where  E  is  the  matrix  of  ones.  Hence  condition  (6a)  follows.  Conversely, 
if  condition  (6a)  is  satisfied  for  some  matrix  E,  then 

n 

IJm(u)  -  m(v)H^  <  max  (  £  E.  .  )||u  -  v|j  . 

l£i£nj«llJ 

Hence  m  is  Lipschltz  continuous.  This  proves  the  proposition. 

The  next  proposition  gives  a  necessary  and  sufficient  condition  for 
(6b)  to  hold.  Recall  that  a  Z -matrix  is  a  square  matrix  with  non-positive 
off-diagonal  entries  and  *  Minkowski  matrix  is  a  Z -matrix  which  is  also 


-13- 


a  P-matrix.  See  Plemmona  [12]  for  various  characterizations  of  a  Minkowski 
matrix.  (Minkowski  matrices  are  also  called  M-matrlces  or  K-matrices  in 
the  literature.) 

Proposition  3.4.  Let  B  and  E  be  respectively,  a  Minkowski  matrix  *nA  a 

--1-  _  _ 
nonnegative  matrix.  Then  p(B  E)  <  1  if  and  only  if  the  matrix  Q  «  B  -  E 

is  a  Minkowski  matrix. 

Proof.  This  result  follows  from  a  characterization  of  a  Minkowski  matrix. 
See  Plemmons  [12]  and  Varga  [14]. 

As  a  consequence  of  Proposition  3.4,  the  next  result  shows  that  a 
necessary  condition  for  (6b)  to  hold  is  that  the  matrix  A  itself  is  an  H- 
matrix  with  positive  diagonal  entires. 

Corollary  3.5.  Suppose  that  condition  (6b)  holds.  Then  A  is  an  H-matrix 
with  positive  diagonals. 

Proof.  Let  A  be  the  comparison  matrix  of  A.  It  is  easy  to  deduce  that 
A  >  B  -  max (DE  ,  jcj)  ■  Q. 

Since  A  is  a  Z -matrix,  the  fact  that  Q  is  a  Minkowski  matrix  ijnplies  that 
A  is  an  H-matrix  (see  Fiedler  and  Pt£k  [6]).  Since  Q  is  Minkowski,  we  have, 
for  each  index  i, 

0  <  «u  <  5lt  -  |cu|  S  Bu  -  clt  -  Au 

because  B  has  positive  diagonal  entries.  This  completes  the  proof  of  the 
Corollary. 

In  the  rest  of  the  paper,  we  derive  several  specialisations  of 
Theorem  3.2.  The  next  result  La  concerned  with  the  nonlinear  complamantarity 


problem  (NLCP) 


(9)  y  £  0  ,  g(y)  >0  and  yTg(y)  -  0. 

Corollary  3.6.  If  by  lotting  m(y)  -  y  -  g(y) ,  there  exiata  a  matrix  E 

with  p(E)  <  1  such  that  condition  (6a)  ia  satisfied,  then  for  an  arbitrary 
0  r  k 

initial  vector  u  ,  the  sequence  [u  ]  defined  by 

(10)  u^1  *  max{0,  m(u^) }  for  k  £  0 
converges  to  a  solution  of  the  NLCP (9). 

Proof.  This  follows  from  Theorem  3.2  by  noting  that  the  iteration  (10)  ia 
a  special  case  of  (1)  under  the  splitting  (I,  0)  of  the  identity  matrix. 

We  mention  in  the  last  section  that  various  iterative  methods  for  the 
LCP(M,  q)  can  be  obtained  from  the  Iterative  scheme  (1)  specialized  to 
the  ICP  (kjM,  k^q,  m)  with  m(x)  »  k^x  under  the  splitting  (B,  C)  of  the 
matrix  kjM  given  In  (4).  The  next  result  is  concerned  with  this 
particular  ICP. 

Corollary  3.7.  Let  and  k2  be  scalars  with  k^  >  0  and  k2  <  1.  Let  D 
and  K  be  matrices  such  that  the  matrix  B  given  in  (4)  is  an  H-matrlx  with 
positive  diagonals.  Suppose  that 

(11)  P(B  1  maxdkjD  ,  JcJ))  <  l 

0  1c 

where  C  ia  given  in  (4).  Then  for  any  initial  vector  u  ,  the  sequence  [u  ) 
defined  by  (1)  specialised  to  the  ICP  (kjM,  k2q,  m)  where  m(x)  -  kjX, 
converges  to  a  solution  of  the  LCP  (M,  q). 

Proof.  This  corollary  is  a  direct  specialisation  of  Theorem  3.2. 
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Corollary  3.6  extends  a  result  of  Aba  [2]  which  deals  mainly  with 
Maagasariaa ' s  Algorithm  II.  Aha  proved  the  result  for  strictly  triangular 
K,  diagoaal  D  sad  kj  ■  0.  The  corollary  also  extends  various  results  of 
Aganaglc  [1].  la  fact,  Agaaagic  required  a  somewhat  stronger  assumption 
than  the  one  (11)  needed  la  the  corollary. 

Specialising  Corollary  3.7  to  the  point  version  of  Algorithm  II,  we 
obtain  the  following  convergence  result  for  a  point  SOR  method  for  solving 
a  LCP. 

Corollary  3.8.  Let  M  be  an  H-matrix  with  positive  diagonal  entries. 

Then,  provided  chat  0  <  <s*  <  2/fl  +  p(Jj)]  where  •  D*1  C|l*|  -  |u|  ]  and 

D,  L  and  U  are  the  diagonal,  strictly  lower  and  upper  triangular  parts  of 

lc 

M  respectively,  the  sequence  {z  }  generated  by  the  point  version  of 
Algorithm  II  converges  to  a  solution  of  the  LCP  (M,  q) .  If  M  is  symmetric 
as  well,  the  same  conclusion  holds  if  0  <  to*  <  2. 

Proof.  According  to  Proposition  2.1  and  Corollary  3.7,  it  suffices  to 
verify  that 

(12)  p[(D  -  <ok|L|)'1(«*|u|  +  |l  -  <s*|D)  <  1. 

Suppose  that  0  <  <u*  <  1.  Then 

(D  -«*|l|)  -  (us* | U |  +  [l  -  u»*|d)  -  <u*(D  -  |L|  -  |u|)  -  «*M 

where  M  is  the  comparison  matrix  of  M.  Hence  by  assumption  and 
Proposition  3.4,  (12)  follows.  On  the  other  hand,  if  a*  >  l,  then  the 
left-hand  term  in  (12)  is  equal  to 

p[(V  -  «*|L|)"1  ((1  -  «*)  D  -  «*fu|)]. 
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By  a  charactarixation  of  H-matrlcaa  (saa  Varga  [13 1)  Che  latter  spactral 
radlua  la  strictly  laaa  than  1  if  <n»  <  2/[l  +  p^)]  whs  re  j  -  d“1[|l|  -  |u[  J. 
Iha  laat  conclusion  of  tha  corollary  la  obvious. 

Va  conclude  thia  papar  by  presenting  a  convergence  raault  for  a  block 
Gauaa-Sldal  method  for  solving  tha  LCP. 

Corollary  3.9.  Lac  M  ba  an  H -matrix  with  poaitiva  diagonal  entrlaa.  Than 

ft. 

tha  aaquanca  {a  }  ganaratad  by  Algorithm  II  with  a*  -  l  convargaa  to  a 
aolutloo  of  tha  LCP  (N,  q). 

Proof*  Lat  0,  L  and  U  ba  aa  apadflad  in  Propoaitian  2.1.  By  Proposition 
2.1  and  Corollary  3.7,  it  suffices  to  show  that  p[(5  -  Jl|)_1|u|]  <  1 
whara  D  is  tha  comparison  matrix  of  D.  Sinca  5  -  [l|  -  |u(  »  M  share  M 
is  tha  comparison  matrix  of  M,  tha  daairad  conclusion  now  follows  from 
Proposition  3.4  and  tha  assumption. 
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In  Part  I  of  this  study,  ve  have  defined  the  implicit  complementarity 
problem  and  investigated  its  existence  and- uniqueness  of  solution.  -^In  the 
present  paper r  me  establishes  convergence  theory^ for  a  certain  iterative 
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nethods\for  solving  a  linear  complementarity  problems - 
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